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Abstract 

The Student's-t mixture model, which is heavily tailed and more robust than the Gaussian mixture model, has recently 
received great attention on image processing. In this paper, we propose a robust non-rigid point set registration algorithm 
using the Student's-t mixture model. Specifically, first, we consider the alignment of two point sets as a probability density 
estimation problem and treat one point set as Student's-t mixture model centroids. Then, we fit the Student's-t mixture 
model centroids to the other point set which is treated as data. Finally, we get the closed-form solutions of registration 
parameters, leading to a computationally efficient registration algorithm. The proposed algorithm is especially effective for 
addressing the non-rigid point set registration problem when significant amounts of noise and outliers are present. 
Moreover, less registration parameters have to be set manually for our algorithm compared to the popular coherent points 
drift (CPD) algorithm. We have compared our algorithm with other state-of-the-art registration algorithms on both 2D and 
3D data with noise and outliers, where our non-rigid registration algorithm showed accurate results and outperformed the 
other algorithms. 
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Introduction 

Point set registration is a key component in many tasks such as 
medical image registration [1], image analysis, computer graphics 
[2] and shape recognition [3]. Iterative closest point (ICP) 
algorithm [4], which assigns correspondences with the closest 
distance criterion, is a popular algorithm for point set registration 
because of its low computation complexity. However, ICP requires 
the initial positions of the two point sets to be adequately close. 
Gold et al [5] proposed the robust point matching (RPM) 
algorithm that performs soft-assignment of correspondences and 
transformation alternately. Chui et al [6] pointed out that the 
processing of alternating soft-assignment of correspondences and 
transformation in the RPM algorithm is equivalent to the 
Expectation Maximization (EM) algorithm with Gaussian mixture 
model (GMM), where one point set is considered as GMM 
centroids and the other one is considered as data. Then Chui et al 
[7] re-parameterized the deformation field using Thin-Plate-Spline 
(TPS) and resulted in the TPS-RPM algorithm. Revow et al [8] 
represented the contour-like point sets using splines and modeled 
them by the probabilistic GMM formulation, where GMM 
centroids were uniformly positioned along the contour. This 
algorithm allows non-rigid transformation between point sets. Tsin 
and Kanade [9] proposed a correlation-based point set registration 
approach, which was later improved by Jian and Vemuri [10] 
(TPS-L2). This approach treats the non-rigid point set registration 
as the minimization of the L2 distance between two Gaussian 
mixture models, and then parameterizes the point sets using 
explicit TPS parameterizations. 



Recently, GMM is widely used to model the non-rigid point set 
registration as it is a natural and simple way to represent the given 
point sets. However, most of the GMM-based non-rigid registra- 
tion algorithms ignore significant effect of the outliers and noise. 
Myronenko and Song [11] introduced the coherent point drift 
(CPD) algorithm, which enforces the points to drift coherently by 
regularizing the deformation field based on the motion coherence 
theory (MCT) [12,13]. The CPD algorithm models the noise and 
outliers using a uniform distribution and uses a manually defined 
parameter co to denote the amount of noise and outliers. However, 
the non-rigid registration is significantly sensitive to CO and it is an 
arduous task to obtain the optimal value of co. Wang Peng [14] 
introduced a hybrid optimization method to estimate the 
parameter co. But, the hybrid optimization algorithm slows the 
non-rigid registration. Christopher et al [15] introduced a 
Student's-t mixture model (SMM) that is a finite mixture model 
based on the Student's-t distribution. The SMM is heavily tailed 
and more robust against noise and outliers than GMM. 
Gerogiannis et al [16,17] proposed a SMM-based rigid registra- 
tion algorithm. However, their algorithms (Christopher et al [15] 
and Gerogiannis et al [16,17]) can only obtain parameters such as 
data mean, covariance matrix and mixing proportion, which limits 
the algorithm to rigid registration. 

It's indicated that, theoretically, GMM is a special case of the 
SMM [15]. Peel and McLachlan [18] modeled point sets with 
outliers and noise for cluster analysis by using SMM. They later 
treated the SMM as a weighted GMM integral form in order to 
obtain closed-form solutions. But unfortunately, they did not 
extend their approach for non-rigid point set registration. To date. 
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it is stUl a challenge to obtain closed-form solutions for the non- 
rigid point set registration by using the SMM [18—20]. In this 
paper, we propose a novel SMM-based non-rigid point set 
registration algorithm. Specifically, we formulate most of the 
registration parameters, such as the deformation field, the equal 
isotropic covariance, and the degree of freedom of the Student's-t 
distribution, and find their closed-form solutions using the EM 
algorithm. The proposed algorithm is computationally efficient 
and more robust against significant amount of outliers and noise 
than existing algorithms. 

The rest of this paper is organized as follows. We introduce the 
main idea of the Student's-t mixture model for non-rigid point set 
registration and details of the proposed method in section 2. 
Section 3 demonstrates experimental results on 2D contour-like 
point sets, 3D surface-like point sets and 3D cloud-like point sets. 
Finally, we present discussion and conclusions in section 4. 



■ D + Ym 



(5) 



where T{*) is a Gamma function, ym is the degree of freedom of 
Student's-distribution for ym, which can change its distribution 
model to fit the data points, and d(x,y,2) = (x-y)T2— l(x-y) is a 
Mahalanobis square distance. We assume Z = a^I, so as to simplify 
the Student's-t distribution for all SMM components. Then the 
SMM takes the form 



(6) 



Methods 

1 Studenfs-t mixture model 

The first point set X^^j)= {xi,...xjy)'^ is the Z)-dimension data 
considered as a data point set. The other point set 
^AfxD^ tVivJVAf)^ is the Z)-dimension SMM centroid set. jVand 
M are the numbers of points in X and F, respectively. The GMM 
probability density function is defined as 



where /7v(a:„;j^,Z) = 



exp ( - ^ {Xn -ymf^ ^ {Xn "J^ 



(1) 



, S is a 



D 1 

covariance matrix, and is a prior probability fory^. satisfies 
the following constraint 



Figure 1 shows the probability density of a univariate Student's t- 
distribution for various degrees of freedom. 

Mathematically, the Student's-t distribution is equivalent to the 
Gaussian distribution when y^°°, and it is also equivalent to 
standard Cauchy distribution when y = 1 [21]. The SMM is more 
robust than the GMM due to its heavy tail for significant amount 
of outliers and noise. According to equation(4), the SMM can be 
represented as a weighted Gaussian mixture model with a Gamma 
distribution and closed-form solutions can then be obtained by 
using the EM algorithm. 

2 Likelihood function for Student's-t mixture model 

We introduce the completion data X^. in EM framework to 
obtain the closed-form solutions. is defined as {X, Z\, 
Z^^ui, Ujy}, where Zmn = {Zn)m- Zmn = 1 if ^« IS Corresponding to 
y^n, otherwise Zmn-^- {i=l,2,...JV) is a random variable 
following Gamma distribution to scale the Gaussian distribution. 
As denoted in [22], we have 



M 
m= 1 



(2) 



M«L„ = i~/rl 



1 1 



(7) 



If the point sets contain outliers or noise, we model both the data 
set and the outliers (or noise) set using two Gaussian mixture 
models as follows 



where 8 is a small value and ^: is a large value. sfy(xn',yfyi,cll) 
represents the outliers and noise in point sets. In brief, we model 
the outliers and noise using the GMM with a different covariance 
matrix. 

A random variable u following the Gamma distribution is 
defined if there are outliers and noise in point sets. Specifically, we 
assume z/~^(a, P), where a and p are set as (X = P = y/2. 
Combining the definition to equation (3), equation (1) can be 
rewritten as 



f(Xn;ym^^m,^) = j Wy^NiXn^ym.^ /u)dH{u) 



(4) 



where H[u) is / distribution. The probability density function of 
the Student's-distribution mixture model is 



\un,zmn = 1 ' '^'^ ' ^ 



(8) 



The likelihood function L by parameterization is denoted as 
follows 



N M 

L{xn]y^,Wm,(y^,y^) =Y^T^Y1 ^mft{xn]yry,,G^,ym) (9) 



The parameter set for non-rigid point set registration is defined 
as ^= (^1, wm^ yivj JVivj JVmj cr^)- We separate 
parameters from each other in the SMM likelihood function by 

L{^')=L{y^,w,„„,<j\y,„)=L{w„) + L{y,„) + L{y,„,G') (10) 

and get the following formulas by combining equations(5), (7) and 
(8) 
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X 



Figure 1. A univariate Student's-t distribution (//=0, <t= 1) for various degrees of freedom. The Student's-t distribution corresponds to a 
Gaussian distribution when y^oo, 
doi:1 0.1 371 /journal.pone.0091 381 .gOOl 



(11) 



N M 



--lm=l 



J2J2^^A- ln(27r) - In 



2^2 



(13) 



3 Estimation of parameters 

We maximize the conditional expectation of parameters in the 
EM algorithm for estimation of parameters of non-rigid registra- 
tion. The conditional expectation in the k+l iteration is given by 

n=\ m=\ 

ln(w(r')/(x„;3^(r'),(a^)<^+",y(r'))) 

where p {y^^^^^ ]Xy^ is a posterior probability density function that 
is denoted as 



(/c+lj _ \ /_ 

m = \ ^ 



(15) 



By maximizing equation (14), we can obtain the update equations 
of M;L'+",i'+",J'L'+"and [u^ f^'K In equation (14), e(^<*+'); 

^*'^'') is a function with respect to 1) (w^; X;^,z^„ = l)and 

E^(^j^ ^ (In Ufn', Xn,Zfnn = ^)- So wc have to calculate 
E^{k+i){um;Xn,Zfnn = ^) and 1) (In w^; x„,z^„ = 1 ) for the 

maximization ofef^(^+^); . Because u is the conjugate prior 



distribution with respect to fy, we can formulate £'^^,(^+1) 

(Ufn', Xn,Zfnn = 1) aS foUowS 



{k) 



(16) 



Using equation (16), £'vj/(^+i) (In w^; x„,z^;^ = 1) is denoted ; 

E^{k+\) (Ig Um] Xn^Z^n = 1 ) = In + ( ^^^^ ) 



In 



yt^+D 



y)n'+D\ ym+D 
where i/^(») is a Digamma function and ^ I I — In 

is a correction for Inu^. 

We separate y, and from each other for equation (14), 
which can then be rewritten as 
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where Qj, Q^, and Q3 are respectively 



(18) 



(19) 



(20) 



M (k+i) (k) ( N 
/^wn /m 



Pmn 



- D ln(27r) - In (cr^^) + i) In 



(A:+l)|| ||2 



In equation (19), the update equation of w^^^^ is 



Jk+l)_ n^\ 



{k+\) 



N 



(22) 



In equation (20), y^^^^ is a solution to 



Jk) 



TV 

E 

n = \ 



2^ /^m« (23) 



Implementation details of equation (21) and how to force the 
SMM centroids move coherently as a group are described in the 
following subsection. 

4 Non-rigid registration with SMM 

The deformation field T of the non-rigid point set registration is 
assumed as T{Y+v) — Y+v{Y)^ where v is the deformation vector 
set. The regularization (p{Y) regularizes the norm of z;, leading to a 
smoothness deformation field that enforces the SMM centroids to 
move coherently as a group to preserve the topological structure of 



the point set [23,24]. Adding the regularization (p{Y) 
to23(y^"^^M^^)^^^^^5^^^^)' equation (21) can be rewritten as 



e3=e3(j,^';^) + ;.^(v) 



(24) 



where X is the weight of the regularization representing the trade- 
off between the goodness of non-rigid registration and regulari- 
zation. After some manipulations similar to [11], Ss^an be 
rewritten as 

a*'"" = - E (-oin(2.)- 1„(.^-)«) 



Pmn 

n=\m=\ ^ 



(<t2) 



(25) 



where GMxM is a Gaussian kernel matrix with elements 

= exp ^— ^ l^^^-^l ^, which reduces the oscillating energy 

at high frequency [24]. G(m,«) is a column vector in the kernel 
GMxM. The parameter (3 defines the width of smoothing 
Gaussian filter in [25]. WMxD is the weight matrix of GMxM. 



(21) Using =0, W^**' can be given by 



W'^'^ = {diag{p'^'^'H)G+X{a^fl) ' 

(p('^+"A:-rfjfl^(p('^+"i)yW) 

where P^^"*"^) is a MxN matrix with elements pi 



(26) 



(^+1) (^+1) 



representing the posterior probability density 

(JA . 

correction by Umn . 1 is a column vector of all ones. I is an identity 
matrix, diag(-) represents a diagonal matrix. Using— ^^-^ = 0, 



,2\(^+l) 



is denoted as 



N M M 11? 



L:„-j<?-G(w,.)WW 



n= 1 m= 1 



(27) 



Theoretically, Xn is considered as an observation in Bayesian 
approach. The discrete latent variable Zmn specifies which 
component of the Student' s-t mixture model generates and 
the continuous latent variable u^^ specifying the scaling for the 
corresponding equivalent Gaussian [15]. Consequently, PFand 
in our algorithm respectively correspond to the ones of the CPD 
algorithm if Zmn — 1 ^^d u^^ — 1 • 

Hence, the main steps in this proposed algorithm are 
summarized as follows. 

Step 1: 



PLOS ONE I www.plosone.org 



4 



March 2014 | Volume 9 | Issue 3 | e91381 



Robust Non-Rigid Point Set Registration 



Initialize the parameter set (y^, w^, cr^,z;), f^, X and the 

convergence threshold t. The initial values of w^'' and (cr^)*^^''are 

N M 



JO), 



:1/M and (a^) 



,2\(0). 



respectively. 



Step 2: 

Construct G: g/,- = exp I 









yf-yf 


- 


21! 



Step 3: 

Calculate the probability density for point sets using equation 
(5). 

Step 3: E-Step. 

Calculate posterior probability density Pmn'^^ using equation 
(15) and Umn^^ using equation (16), respectively. 
Step 4: M-Step. 

Estimate Wm^^^ using equation (22), y'm^'^^ using equation (23), 
W^'^'^^^ using equation (26), and (cr^)^^^^^ using equation (27), 
respectively. 

Step 5: 

Calculate the new positions of the SMM centroids using 
Step 6: 

Return to step 3 if the registration error e>x. 

Results 

To compare our SMM-based algorithm with CPD and TPS- 
RPM for non-rigid point set registration, we performed experi- 
ments on 2D contour-like point sets with outliers and 3D surface 
point sets with noise. 



1 Qualitative evaluation 

Here we present qualitative comparisons of the non-rigid 
registration algorithms. The experiments were conducted on 2D 
contour-like point sets, which were from the boundary of 
cerebrospinal fluid (CSF) segmented from CT images. 

In the first experiment, we segmented CSFs very carefully and 
got a quite ideal segmentation result. The corresponding point sets 
are shown in Figure 2(a). The black point set had 488 points and 
the red one had 398 points. The two point sets were considered as 
the ideal boundary of CSF without any outliers or noise. However, 
the boundary in the central region (as indicated by the arrow in 
Figure 2(a)) was tortuous, which caused the CPD algorithm to fail 
to align correspondences between two point sets. Consequently, 
the deformation field obtained by the CPD algorithm was aliasing 
and the topological structure was changed. In Figure 2(d), several 
correspondences in the margin are improperly aligned by using 
TPS-L2. Most of red points were aligned accurately by the TPS- 
RPM algorithm due to its TPS regularization. But because of its 
explicit TPS regularization, the boundary of the central region is 
too smooth to fit data points. Our SMM-based algorithm aligned 
the central points more accurately than TPS-RPM and the 
deformation field produced by our algorithm was smoother than 
the one produced by the CPD algorithm. We set P — 2.2 and A = 2 
in the qualitative experiments. The deformation vectors (as 
indicated by the green vectors) are also displayed in Figure 2(b) 
~ Figure 2(e). 

The deformation vectors produced by the four algorithms are 
zoomed in for the central points in Figure 3. All points were 
aligned properly by our algorithm (see Figure 3(a)). However, 
some points were aligned improperly (as indicated by the arrows a 
and b in Figure 3(b)) by the CPD algorithm, and several vectors 
were crossed (as indicated by the arrows c and d in Figure 3(b)), 
which broke the topological structure of the aligned point set. The 
deformation vectors produced by TPS-L2 are more regularized 





(a) 



(b) 




Figure 2. Registration of 2D CSF boundary points without outliers or noise. The 2D CSF points are segmented on the images from http:// 
www.insight-journal.org/rire/download.php. (a) Ideal CSF boundary points, the black point set contains 488 points and the red one contains 398 
points, (b) The registration of our SMM-based algorithm {p = 22, X = 2), (c) CPD, (d) the Jian and Verimu's algorithm [10] (TPS-L2), (e) TPS-RPM. The 
green vectors denote the deformation vectors produced by non-rigid registration algorithms. 
doi:1 0.1 371 /journal.pone.0091 381 .g002 
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(c) (d) 

Figure 3. Deformation vectors of central points, (a) Our 

algorithm, (b) CPD, (c) TPS-L2, (d) TPS-RPM. The deformation vectors 
produced by our algorithm are more regularized than the ones 
produced by CPD, and the registration result of our algorithm is more 
accurate than TPS-L2 and TPS-RPM. Some deformation vectors 
produced by CPD were crossed, which broke the topological structure 
of the point sets. The aligned point set was too smooth to fit the CSF 
boundary points by using TPS-RPM. 
doi:1 0.1 371/journal.pone.0091 381 .gOOB 



than the ones produced by our SMM-based algorithm and the 
CPD algorithm. However, the TPS-L2 algorithm failed to align a 
few correspondences in the central region (as indicated by the 
arrow in Figure 3(c)). Similar to TPS-L2, the deformation vectors 
produced by the TPS-RMP algorithm are regularized, but TPS- 
RMP also failed to align the central points (as indicated by the 
arrow in Figure 3(d)). Obviously, our SMM-based registration 
algorithm is the most accurate one of the four algorithms. 

Generally, segmentation results are sensitive to image quality 
and the selection of parameters of segmentation algorithms. 
Consequently, we may not segment targets very accurately and 
obtain ideal point sets without outliers and noise. In the second 
experiment, the comparison of the registration algorithms was 
performed on the point sets of CSF boundary containing outliers 
(see Figure 4(a)). There were 184 outliers clustered into 1 1 outlier 
sets in the red point set and 54 outliers clustered into 4 outlier sets 
in the black point set. The registration result by our algorithm is 
shown in Figure 4(b). Figure 4(c) (the registration result by the 
CPD algorithm) shows that one outlier set (as indicated by the 
arrow) are shifted towards the CSF boundary, which implies the 
CPD algorithm cannot align the correspondences between the 
point sets when the counter-like point sets contain outliers. 
Figure 4(d) shows the registration result (which is the worst) by the 
TPS-L2 algorithm, indicating that outliers could seriously affect 
the TPS-L2 algorithm. Similar to Figure 4(d), Figure 4(e) also 
shows a bad registration result by the TPS-RPM algorithm, 
indicating that the TPS-RPM algorithm failed to parameterize the 
deformation field of CSF boundary exactly. 

Figure 5 shows the iterations of the four non-rigid registration 
algorithms. The first, second, third and fourth rows demonstrate 
the iterations of our SMM-based algorithm, the CPD algorithm, 
the TPS-L2 algorithm and the TPS-RPM algorithm, respectively. 





(c) (d) (e) 

Figure 4. Registration of CSF boundary point sets with outliers, (a) CSF boundary point sets with outliers, 184 red outliers clustered into 1 1 
outlier sets, and 54 black points clustered into 4 outlier sets, (b) Our algorithm, (c) CPD, (d) Jian and Verimu's algorithm [7] (TPS-L2), (e) TPS-RPM. One 
outlier set was aligned to the CSF boundary points by CPD, and two outlier sets are aligned to the boundary points by TPS-L2 and TPS-RPM. 
doi:1 0.1 371 /journal.pone.0091 381 .g004 
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The degree of the freedom y can change the Student' s-t 
distribution model to fit the data points [26], making the 
convergence of our algorithm only took 40 iterations, less than 
the CPD algorithm, and significantly less than the TPS-L2 
algorithm and TPS-RPM algorithm. 

2 Quantitative evaluation 

2.1 3D face point set registration. In this experiment, we 
demonstrate the quantitative performance of the algorithms on 3D 
face data with 392 points (which are from http://www.csee.ogi. 
edu/myron/matlab/cpd/). We define the registration error as 



E \\^n-yn\\a, 



n=\ 



ifter 



N 

E \\^n-yn\ 
n = l 



X 100% 



(28) 



before 



where ||^« — J«||a/?er denotes the Euclidean distance of the 
correspondences in X and Y after registration, \\Xn—yn\\before 
denotes the Euclidean distance of the correspondences before 
registration, and N is the amount of points in each point set. 

Figure 6 shows the non-rigid registration results produced by 
our SMM-based algorithm, CPD, TPS-L2 and TPS-RPM on the 
3D face data, respectively. In the experiment, we set ^ — 2.2 and 



A = 3, respectively. The registration error is 1.27% by using our 
SMM-based algorithm, and the errors are 2.53%, 3.11% and 
3.54% by using CPD, TPS-L2 and TPS-RPM respectively. 

Figure 7 shows the registration results on the 3D face data 
containing 40% artificially added Gaussian noise with cr = 5. It is 
challenging for non-rigid registration algorithms to align 3D points 
with such noise. Figure 8 shows the registration results on the 3D 
face data, hiding the artificially added noise. Figure 8(a) demon- 
strates that our SMM-based algorithm can align the correspon- 
dences precisely. Figure 8(b) and Figure 8(c) show that CPD and 
TPS-L2 are able to align the point sets accurately except few 
points in the margin (as indicated by the arrow in Figure 8(b) and 
Figure 8(c)). Figure 8(d) shows a registration with serious error by 
TPS-RPM, indicating a failure of TPS-RPM for handling the 
marginal points with so much noise. In Figure 7 and Figure 8, the 
registration error is 3.36% by using our algorithm, and the errors 
are 10.7%, 15.4% and 21.1% by using CPD, TPS-L2 and TPS- 
RPM respectively. 

A series of experiments demonstrate that the registration of the 
CPD algorithm is sensitive to the parameter w which indicates the 
amount of noise in the point sets. However, this parameter had to 
be set manually. We tested the impacts of co on 3D face data with 
40% Gaussian noise and 40% uniform noise, respectively. Figure 9 
shows the great impacts of co on registration. In practice, the 
relationship between the optimal value of (D and the amount of 



SMM- 
based 





CPD 





TPS- ' 
L2 





TPS- 
RPM 



















Iteration 10 



Iteration 20 



Iteration 30 



Iteration 40 



Iteration 50 



Figure 5. Iterations of the four non-rigid registration algorithms. The convergence of our algorithnn takes 40 iterations, while the CPD 
algorithnn takes about 50 iterations and the TPS-L2 algorithnn and the TPS-RPM algorithnn takes nnore than 50 iterations. 
doi:1 0.1 371/journal.pone.0091 381 .g005 
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■■'V (M.'^' 



(a) 



(b) 



(c) 



(d) 



(e) 



Figure 6. Registration of 3D face data sets with 392 data points, (a) 3D face data set, (b) our algorithm {p = 22, 1 = 3, e,= 1.27%), (c) 
CPD(e, = 2.53%), (d) TPS-L2(e, = 3.1 1%), (d) TPS-RPIVl(e, = 3.54%). 
doi:1 0.1 371 /journal.pone.0091 381 .g006 

noise is nonlinear. It's thus a challenge to get an optimal value of respectively. The results demonstrate that our algorithm is more 
CO. Figure 10 shows the registration errors of the four algorithms on accurate and robust than the other algorithms, especially in 
the 3D face point sets with Gaussian noise and uniform noise, presence of significant amount of noise. 





(a) 



(b) 






(c) 



(d) 



(e) 



Figure 7. Registration results on 3D face data containing 40% artificially added Gaussian noise with 6=5, (a) 3D face data sets with 
Gaussian noise {(t = 5), (b) our SMM-based algorithm(e, = 3.36%), (c) CPD(e,= 10.7%), (d) TPS-L2(e,= 15.4%), (e) TPS-RPIVl(e,= 21.1%). 
doi:1 0.1 371/journal.pone.0091 381 .g007 
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(a) 



/ • °. • o ^ • 



• • i 



5r- °. 'o •> 



(b) 




on 5 different 3D lung datasets from 5 different subjects (which are 
from http://www.dir-lab.com/index.html). Each dataset has a 
pair of 3D lung point sets, one set was extracted from the 
maximum inhalation phase image and the other set was extracted 
from the maximum exhalation phase image. Each 3D lung point 
set has 300 cloud-like points, which were selected by experts to 
make the two point sets correspond to each other. It is a herculean 
task for non-rigid registration algorithms to match cloud-like point 
sets accurately due to lack of the topological structure or the 
geometry structure in such sets. For the above reason, we did not 
add any noise or outliers in these cloud-like point sets. The average 
registration errors are 0.26 mm, 1.05 mm, 0.66 mm and 2.97 mm 
by using our SMM-based algorithm, CPD, TPS-L2 and TPS- 
RPM respectively. Figure 1 1 shows the non-rigid registration 
results produced by our SMM-based algorithm and the other 
three algorithms. Figure 1 1 (b) shows that our algorithm aligned all 
correspondences accurately. Figure 1 1 (c) and Figure 1 1 (d) show 
that CPD and TPS-L2 are able to align most correspondences 
except few points in the bottom of the lung. Figure 1 1 (e) shows a 
bad registration result by using TPS-RPM, indicating a failure of 
TPS-RPM for handling these cloud-like points. Figure 12 shows 
the registration errors of the four algorithms on the five subjects. 

Discussion 



(c) 



(d) 



Figure 8. Registration results on 3D face data without the 
display of artificially added noise, (a) Our SMM-based algorithm (b) 
CPD, (c) TPS-L2, (d) TPS-RPM. The correspondences are aligned 
accurately by our algorithm, which demonstrates that our SMM-based 
algorithm is robust against the significant amount of noise. The CPD 
algorithm and the TPS-L2 algorithm failed to align some correspon- 
dences in the margin, and the TPS-RPM algorithm absolutely failed to 
align the marginal face points. 
doi:1 0.1 371 /journal.pone.0091 381 .g008 

2.2 3D lung point set registration. In this experiment, we 
demonstrate the quantitative performance of the four algorithms 



The Student' s-t mixture model is heavily tailed and more robust 
than the Gaussian model, and has recently received great attention 
on image processing. We have proposed a probabilistic method for 
non-rigid point set registration based on the Student' s-t mixture 
model. We considered the alignment of two point sets as 
probability density estimation, where one point set was represent- 
ed as the Students'-t mixture model centroids, and the other one 
was represented as the data points. We iteratively aligned the 
SMM centroids by the EM algorithm and estimated the posterior 
probabilities of centroids, which provided the correspondence 
probability. Finally we estimated all registration parameters and 
obtained closed-form solutions. It's worth noting that, in the 
Student's-t mixture model, we formulate all registration param- 
eters of each point using the expectation maximization (EM) 
algorithm, making the coherent point drift algorithm modeled by 
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Figure 9. Impacts of o) on registration of CPD. The registration results of CPD are sensitive to the parameter co. 
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Figure 10. Registration errors on the face data with noise, (a) Registration errors on face data with Gaussian noise, (b) registration errors on 
face data with uniform noise. Our algorithm is robust against uniform noise and Gaussian noise, and outperformed CPD and TPS-RPM. 
doi:1 0.1 371/journal.pone.0091 381 .gOI 0 



the Gaussian mixture model be a special case of our algorithm. 
Our algorithm is superior to the other algorithms for addressing 
the non-rigid point set registration problem, especially when 
significant amounts of noise and outliers are present. 

Our contribution includes the following aspects. Firstly, we 
modeled the non-rigid point set registration using the Student's-t 
mixture model (SMM). Our SMM-based algorithm is more robust 
than CPD, TPS-L2 and TPS-RPM when the point sets contain 
significant amounts of noise and outliers. Secondly, we formulated 
all registration parameters of each point using the EM algorithm, 
making the CPD algorithm to be a special case of our algorithm. 
Thirdly, we estimated the prior probability Wj^ in the EM 
algorithm while it was assumed as a constant in the CPD 



algorithm, which makes our algorithm robust against outliers and 
noise. Finally, we used the EM algorithm to estimate all 
parameters and obtained closed-form solutions except y, making 
our registration algorithm computationally efficient. 

The data used in preparation of this article are available to the 
public. The 2D CSFs point sets are segmented from the images 
from the Retrospective Image Registration Evaluation Project 
(http://www.insight-journal.org/rire/index.php). The 3D face 
data can be loaded from http://www.csee.ogi.edu/myron/ 
matlab/cpd/. The 3D lung point sets can be obtained from the 
Deformable Image Registration Laboratory (http://www.dir-lab. 
com/index.html.). 
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Figure 1 1. 3D lung point set registration on a subject, (a) 3D lung point sets before registration, (b) SMM-based algorithm, (c) CPD, (d) TPS-L2, 
(e) TPS-RPM. Our algorithm performs the best. 
doi:1 0.1 371 /journal.pone.0091 381 .gO1 1 
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Figure 12. Performance of the four non-rigid registration algorithms on 5 subjects. Our algorithm is the most accurate one of the all 
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